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Abstract. A very recent observation by the Auger Observatory group [1] claims 
strong evidence for cosmic rays above 5.6 xlO 19 eV - (56 EeV) being protons from 
Active Galactic Nuclei. If, as would be expected, the particles above the ankle at about 
2 EeV are almost all of extragalactic (EG) origin then it follows that the characteristics 
of the nuclear interactions of such particles would need to be very different from 
conventional expectation - a result that follows from the measured positions of 'shower 
maximum' in the Auger' work. 

Our own analysis gives a different result, viz that the detected particles are still 
'massive' specifically with a mean value of (In A) = 2.2 ±0.8. The need for a dramatic 
change in the nuclear physics disappears. 



PACS numbers: 96.50.sb, 96.50.sd, 98.70.Sa 
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1. Introduction 




Figure 1. AUGER source map showing possible 'sources' A, B and C [1]. The energy 
threshold is 57 EeV. 

It has long been suggested that the particles above the ankle are extragalactic (e.g. 
[2, 3]); indeed, some believe that the transition starts at an even lower energy than 2 
EeV (e.g. [4]). There have been many claims for EG 'signals' from specific sources (e.g. 
[5]) but, apart for rather strong evidence for particles from the VIRGO cluster (the 
centre of the supercluster in which we are situated) the results have been conflicting. 
There were thus high expectations for the results from the very large Auger Observatory 
and such result, based on an exposure (area times time) exceeding the sum total of the 
world's data have recently appeared [1]. Figure 1 shows the results and it is evident that 
there is 'dumpiness' in the arrival directions. The authors draw attention to a number 
of features, principally 

1 The presence of a cluster of events round the CEN-A radio source. 

2 Coincidences, above the chance level, with known Active Galactic Nuclei (AGN) 
out to 75 Mpc. 

Their conclusion that the primaries are protons is based on the contention that 
the deflections in the intergalactic medium (IGM) and the Galaxy would nullify the 
coincidences. 

Although not stated, the need for a change in the Nuclear Physics follows from 
examination of the world's data (and their own e.g. [6, 7]) on the depth of shower 
maximum, which indicate (In A) ~ 1.4 at 10 EeV and ~ 2.5 at 40 EeV, the highest 
energy point plotted in the Auger results. With the conventional Nuclear Physics model, 
protons ((In A) = ) are certainly ruled out for the particles above 56 EeV. If true, 



Heavy Cosmic Ray Nuclei from Extragalactic Sources above "The Ankle" 3 

this result would arguably be more important than the demonstration that AGN are 
responsible for the ultra high energy particles. 

This, then is the problem addressed here: Are 'heavy nuclei' ((In A) — 2 3) ruled 

out? 

2. The Ankle 

As remarked already, and referred to by us in several publications (eg. [8]) we consider 
that this feature marks the transition from a mainly Galactic (G -) to a mainly 
Extragalactic (EG -) origin. Some others have it as a property of EG protons and 
a demonstration that this is the case would clearly support the Auger contention. We 
have made many arguments against the EG protons/ankle hypothesis (eg. [8]) and 
these are strengthened from observation of the Auger energy spectrum reported in [7]. 
The ankle is so sharp as to make its explanation in terms of E - p quite untenable. 

The situation can be seen by reference to Figure 2. A two component spectrum with 
the spectra Galactic (G -) and Extragalactic (EG -) crossing sharply, as in Figure 2a, 
clearly gives a sharp ankle. The EG protons alone, spectrum (Figure 2b) has too 
smooth a transition; in fact, various factors would make the transition smoother still for 
such a model [8]. Figure 2c refers to one of our variants [8]; specifically, origin of ultra 
high energy cosmic rays (UHECR) in quasars. It is our contention that the actual form 
of the spectrum can be used to define the best fit spatial distribution of the sources as a 
function or red shift Q(z) see [8]. The rather sharp fall in the Auger spectrum at about 
40 Eev has relevance here, as will be discussed in more detail elsewhere. At present 
a possibility is that heavy nuclei are involved, indeed,this is the basis of our following 
arguments. In fact the prediction shown,which is for protons alone would need to be 
displaced downwards to allow for 'heavy' nuclei at energies below 20 EeV,too. 

3. The Depth of Shower Maximum 

Measurements over many years have shown that the depth of maximum increases with 
increasing energy and its value is roughly mid way between expectation for 'all protons' 
and 'all iron' for essentially all the nuclear models to date. With their superior statistics 
and analysis, the Auger work [7] has shown that there is structure in the energy 
dependence, with a feature near the ankle energy. Figure 3 shows the results and 
Figure 4 shows the resulting (In A) from our analysis. It is interesting to note that 
the HiRes EAS array shows a similar feature: an Xmax change close to the ankle [9]. 
A relevant matter to consider now is the expected mass composition before the ankle. 
This cannot be anywhere near mainly 'protonic' because of the lack of large anisotropies 
favouring the Galactic Plane, as shown by us in a previous detailed analysis [10]. Thus, 
the Nuclear Physics models should not be too inaccurate here. Were the particles to 
be mainly protonic only above the ankle, the change in Nuclear Physics Model would 
need to take place over half a decade of energy, at most, and we consider this to be 
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Figure 2. The AUGER energy spectrum [7] in comparison with various predictions 

a. Our model fit [8] where the Galactic (G -) and Extragalactic (EG -) spectra are 
simple power laws. 

b. Comparison with EG protons only model of [4]. 

c. Comparison with our Q2 model for protons [8]. The shortage above 10 EeV would 
be covered by heavier nuclei. 

unphysical. 

At this stage, it can be remarked that, in fact, the Auger Amax results would give 
too high an anisotropy at lEeV,where the particles are of Galactic origin. The lower 
Amax values measured by most others would not [10]. 
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Figure 3. Depth of shower maximum versus energy measured by different experiments 
compared with different model predictions. 



4. The arrival directions 

4.1. CEN-A 

Returning to Figure 1, together with the Auger authors, we are impressed by the signal 
from Centaurus-A (CEN-A), a long favoured source with its double jet, high power and 




E (eV) 

Figure 4. (In A) vs energy from our analysis of the AUGER results and different 
models. Most other X max values from Figure 3 would give higher mean masses. 
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flat radio spectrum. We would argue that the 8 particles within 20% of CEN-A are 
probably due to it. The lack of better symmetry about the source can be understood in 
terms of the varying collection efficiency in this region together with lack of symmetry 
in the magnetic fields in the IGM. We call this region A. Two more regions may be 
populated by specific nearby sources, denoted B and C in the Figure. 

Although it is claimed that there is an excess above chance of coincidences with 
AGN in general the statistics will be made worse when the CEN-A events are removed, 
and, furthermore, there is a remarkable lack of events in the region above CEN-A, 
where the number of AGN is seen to be very large. This approaches the VIRGO region 
(already mentioned), where the cluster is situated 17 Mpc away. It is instructive to make 
an estimate of how many particles might have been expected to be seen by Auger from 
VIRGO. Including the difference in collection efficiency, a factor ~ 2.5, we would expect 
to see, for a single CEN-A source at the distance of VIRGO, about 0.3 events,therefore 
there are less than a few CEN-A type galaxies amongst several thousand galaxies, and 
probably ~ 10 AGN in that cluster. Were all AGN like CEN-A we would expect to 
see ~ 3 events here. This, not unexpected, variability of output in CR amongst AGN 
of different types coupled with the lower detection efficiency of distant AGN and not 
to mention their time variability tends to give further doubt to an analysis in which 
coincidences are sought between a non homogeneous set of AGN and UHECR. In fact, 
there was a prior likelihood of large radio galaxies rather than AGN in general being 
likely sources or UHECR for two reasons: 

a. the likely detection of M87, a radio galaxy with a pronounced jet in the Northern 
hemisphere, [3] and, 

b. the obvious need for a type of source with large linear dimensions and, preferably 
radio spectra with small exponent indicating, for electrons at least, flat energy 
spectra. 

For radio galaxies, those giving the highest radio fluxes at earth are, in order of increasing 
distance from earth: 

• CEN-A (NGC 5128) at 4.9 Mpc 

• VIRGO A (NGC 4486, M87) at 16.8 Mpc 

• and FORNAX (NGC 1316, ARP 154) at 16.9 Mpc 

4.2. 'Source B' 

In the list just given, the first two have been mentioned already. FORNAX is seen 
to be not far from the 'Source B' but probably too far to be physically associated. 
However, it is in the FORNAX cluster and this has galaxies extending across to /, b: 
200°, —40°. Most notable is that of the radio sources with flat radio spectra (associated 
with elliptical galaxies), [13] the flattest, with exponents ~0 and 0.1, are near to Source 
B. They are NGC 1052 and NGC 1407 at I, b, distance: 182°, 58°; 17.8 Mpc and 209°, 
-50°, 21.6 Mpc. 
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We conclude that there are reasonable contenders for 'Source B'. 

4.3. 'Source C 

The evidence for this 'cluster' of UHECR arrival directions being associated with a 
single known source is not strong. There are no obvious candidates. There are only a 
few 'normal' galaxies within 20 Mpc [14, 15] although there is a nearby galaxy within 5 
Mpc. Presumably a source further away is responsible? 

4-4- Other Source Complications 

Its well known that many AGN are time variable (and CEN-A is no exception). Thus, 
in view of transit time differences between UHECR and protons for very distant sources, 
the optical and UHECR sky may differ. It is interesting to note that nearby colliding 
galaxies (some of which go on to produce AGN) have not (yet) been seen (see [15] for 
previous work). In addition, to the different types of AGN their distances are clearly of 
great relevance; thus, catalogues are needed of putative claims for coincidences giving 
particle energies and AGN distances (and types). It can be remarked at this stage 
that the term 'Active Galactive Nuclei' is probably a misnomer to describe 
the UHECR sources. The large radio sources have ceased to have active nuclei 
by the time the radio jets are seen. 

5. The case for, or against, non protons 

5.1. Acceleration Mechanisms 

Starting with acceleration, there is an obvious advantage in accelerating high Z particles 
insofar as the commonly considered acceleration mechanisms operate, with a rapidly 
falling energy spectrum, to some maximum rigidity. Thus, provided there is an adequate 
pool of pre-accelerated nuclei, say up to iron, such non protons should be at an 
advantage. There is a similar situation in the tens of thousands of GeV region, where 
the mean mass increases with energy [16]. 

5.2. Angular deflections in the Galaxy 

In [1] it is claimed that the deflections in the Galactic magnetic field do not allow the 
incoming nuclei to maintain their common direction. There are two points that should 
be made here: 

a. For CEN-A the trajectories are very largely along the coherent magnetic field 
direction (the Sagittarius Arm). 

b. Calculations for the random field, using the parameters of field strength and reversal 
length [17, 18], yield a median deflection of 0.72 degrees, again for CEN-A 
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Figure 5. Integral distribution of angular deflections (6 - in degrees) for protons from 
CEN-A (line labeled 5 Mpc) and the distributions for 15 Mpc and 100 Mpc distant 
sources for a comparison. 



Turning to 'Source B' at I ~ 190°, b ~ 60°, the direction gives paths mainly in our weak 
'spur'. 

Finally, for 'Source C at / ~ 60°, b ~ 45°, again the direction is along a spiral arm, 
specifically along the inter arm region between the Orion and Sagittarius arms where 
the magnetic field is lower, as well as being mainly ineffectual because of direction. 
Furthermore, the length of path through the irregular field is smaller than that for 
CEN-A (by a factor sin sin 19°/ sin 45° = 0.46). 

5.3. Angular deflections in the IGM 

The problem of calculation is more severe here, largely through uncertainty in the 
reversal length of the field, or to be more specific, the power spectrum of the irregularities 
[19]. Elsewhere [15] we gave an approximate expression for the root mean square angular 
deflection 5105 ^ (da) (pc) degrees, where d is the distance to the source in Mpc, a is 
the reversal length in Mpc and pc is the energy in units of 1000 EeV. B is in the field 
in /iG. Here, we have pc ~ 0.7. The result is similar to what we find using our more 
accurate method; the actual integral distribution is given in Figure 5 for the case of 
CEN-A (d = 5 Mpc); here, the median deflection is 1.1° for protons. 

6. Application to the UHECR Map 

Our method is to use the order of magnitude values of the median angular deviations 
from Figure 1, to give the expected median values of Z. Converting the mean value 
for the sources to an effective mass (A = 2Z) gives our estimated (In A). This is then 
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Table 1. Median expected displacements (in degrees) for protons from the sources 
indicated, and their 'total', i.e. addition in quadrature. Comparison with observed 
displacements gives an order of magnitude estimate of the particle charge, Z. 



Source 


Distance (Mpc) 


Galaxy 


IGM 


Total 


Median displacement observed 


Z 


CEN-A 


5 


0.7 


1.1 


1.3 


10 


7.7 


Source B 


20 


0.46 


2.2 


2.2 


6 


2.7 


Source C 


33 


0.48 


2.8 


2.8 


10 


3.6 



indicated on Figure 4. 

The value for CEN-A, the best identified source, is seen to be (Z) = 7.7 and (In A) 
follows as ~ 2.7. Taking the mean of all three gives (Z) = 4.7 and (In A) = 2.2. It seems 
to us unlikely that the true value is outside then limits; certainly, (In A) = appears 
not to be needed. 

7. Conclusions 

We conclude that it is probably not necessary to change the Nuclear Physics of high 
energy interactions at energy above 60 EeV, or so. 

The way forward in the analysis of the Auger results is to endeavour to check the 
hypothesis that 'nearby' (within some 10s of Mpc) flat spectrum radio galaxies are 
responsible. Identification will clearly rely on examination of the allotted energies to 
events within clusters as a function of radial distance from the possible source. Individual 
Amax values need treating in the same manner. 

A complication, affecting all searches, is the fact that the distant source may not 
be seen optically to be 'still on' when the particles arrive, [20]. 
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